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I If Be solar neutrinos can be observed in real time experiments, then an 

(N ■ 

^ ' extended region of Am can be probed by a proper analysis of the rapidly 

a^ : 

I changing phase of vacuum oscillations due to the eccentricity of the earth's 

I orbit about the sun. For the case of maximal vacuum mixing, a kind of Fourier 

' analysis of expected data for one year's time could uniquely pick out Am^ if 

PLiI it lies in the region --^ 10~^° — 6 x 10^^ (eV)^. 

PACS: 14.60.Pq, 13.10.+q, 25.30.Pt 

X ■ I. INTRODUCTION 

In a previous work it was argued that with maximal vacuum mixing there is agreement, 
with minor modifications, between extant observations of solar neutrinos and predictions by 
the standard solar model (SSM) The maximal vacuum mixing case considered was 

that in which the phase of neutrino oscillations coming from the sun is averaged, leading to 
50% of the neutrinos arriving at the earth as electron neutrinos. As a result of this averaging, 
while sin^ 29 was assumed to be maximal (equal to one), Am^ was not determined and taken 
to lie in the approximate range 10~^ < Am^ << 10^^ with an exclusion of the approximate 
range 3 x 10^^ < Am^ < 10~^ for maximal mixing |^ due to the lack of an observed 
day-night effect in the SuperKamiokande data 0. 
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On the other hand, the recent first results of the SNO measurement of charged current 
interactions produced by neutrinos ||^, taken in combination with the elastic scattering 
result of the Super-Kamiokande collaboration [|^, indicate that only about one third of the 
neutrinos arriving at the earth from the sun are electron neutrinos, with the other two thirds 
being /i or r neutrinos. Oscillation into sterile neutrinos now seems relatively unlikely from 
the SNO result. 

While at first glance this comparison seems to make maximal vacuum mixing less likely, 
a global analysis of the SNO result with the other solar neutrino experiments, chlorine 
0, Super-Kamiokande 0, and gallium p!0| , pT| has led to the conclusion that "the CC 
measurement by SNO has not changed qualitatively the globally allowed solution space 
for solar neutrinos, although the CC measurement has provided dramatic and convincing 
evidence for neutrino oscillations and has strenghened the ths case for active oscillations 
with large mixing angles |jl2[-" Furthermore, global analyses [0 [|13| do not completely 



exclude solutions to the solar neutrino problem in the mass region 10"^*' < Am^ < 10^*^ for 
maximal (or near maximal) mixing. In the following, the time varying phase of oscillating 
^Be neutrinos is investigated as a possible method to discover (or exclude) a solution of the 
solar neutrino problem in that mass region. 

In the mass region 4x 10^^^ < Am? < 10^^ there are so-called "just-so" vacuum solutions 
of the solar neutrino problem, where the phase of the oscillation of ®B neutrinos coming from 
the sun is not completely averaged |1^,|15|]. Recall also that there is a large change 



in the ^Be electron neutrino flux over the year in the ^B "just-so" region due to the change 
in phase of order 7r/2 in a year brought about by the ±1.67% yearly orbital variation from 
the mean distance of the sun to the earth. As will be shown in the following, when phase 
averaging due to the temperature of the sun and phase damping due to the MSW effect 
are considered, it turns out that phase variation in ^Be neutrinos should be observable for 
Am? in the range from about 10"^'^ up to about 6 x 10^^ (sV)^. This observable range of 
Am^ via ^Be neutrinos turns out to be in approximate agreement with a previous analysis 



by de Gouvea, Friedland, and Murayama |]18 



II. OSCILLATIONS: THERMAL AVERAGING; MSW DAMPING 



There is a low energy region of solar neutrinos dominated by the nearly monenergetic 
862 KeV line from electron capture by ''Be in the sun. It is the purpose of the Borexino 
19| experiment, soon to come on line, to measure these neutrinos in real time. Feasability 



studies are also being carried out for other experiments to measure ^Be neutrinos in real 



time such as LENS [|0[ and HELLAZ |T| 



If the ''Be neutrinos were truly monoenergetic then the number of neutrinos detected via 
electron scattering in an experiment like Borexino (normalized to unity for no oscillations) 
would take the following form for vacuum oscillations 



Am^) = 1 - 0.79 sin^ 26 sin^ ^^"^'-^W n ) 

^ (0.862)0.00248' ^ ^ 

where 6 is the vacuum mixing angle, Am^ is expressed in (eV)^, the /i or r neutrino scattering 

relative to electron neutrino scattering at 0.862 MeV is 0.21 ||2^, and 



m = X 10«, (2) 

1 + e cos 

the distance from the Earth to the center of the sun (in km.), which varies through the year 
due to the eccentricity e = .0167 of the Earth's orbit about the sun. Note that since we 
take the number of neutrinos detected as a function of 0, the phase of the earth in its orbit 
about the sun, rather than of the time of the year there is no 1/L^ seasonal variation in 
R; it is canceled by the Jacobian in going from time as an independent variable to as an 
independent variable (Kepler's second law). 

It has been pointed out by Pakvasa and Pantaleone that the ^Be energy line is 



thermally broadened not only by the spread in nuclear velocities (Doppler broadening) but 
also by the solar temperature of the approximately 80% of the capture electrons that come 



from the continuum. We makes use of the published table of Bahcall |^ which was obtained 
by convoluting both these sources of thermal spreading to obtain the energy profile of the 862 
keV ^Be solar neutrino shown in Figure 1. Note that the distribution Wt{x) is asymmetric 
in shape and plotted as a function of x = {E — 0.862)/0.862 in MeV. 



Expressing Am^ in 10 ^(e^)^ one obtains 



RU, Am") = 1 - 0.79 sin' 29 [ dx Wt(x) sin' (- 2198Am' ^\ 

^ J ' V l + .0167cos(/) l + x ^ ^ 



0167cos(/))(l + a;);' 



or 

4396Ato' 



Am') = 1 - 0.395 sin' 29 



1 — dx Wt{x) cos I 



(4) 



, (1 + . 0167 cos 0)(l + x) 

For clarity and convenience we will retain these constants explicitly and express Am' in 
units of 10~^(ey)' unless otherwise specified for the rest of this paper. 

Since x is constrained to contribute to Eq.(9) only when it is much smaller than unity, 
one may set + x) equal to 1 — x and carry out the integration to obtain 



Am') = 1 - 0.395 sin' 29 
-I 



^ , Am' \ / 4396Am' . 

1 cos — — — (5) 



'■s 



1 + .0167 cos 07 V 1 + .0167 cos 

Am' \ / 4396Am' 
jsm(^- 



1 + .0167 cos (f)J V 1 + .0167 cos 



with 

In additon to temperature damping of the the oscillations there is a second factor that 
we may call "MSW damping" . With maximal mixing, the phase averaged rate of electron 
neutrinos does not depend on whether an MSW transition has taken place. However, if an 
MSW transition has taken place in the sun, then the maximally mixed neutrino emerges 
in the form of a pure mass eigenstate (i.e. Ai = 0, A2 = 1 in Eq.(A5) of Appendix A). 
Although the pure mass eigenstate > is half electron neutrino and half other flavor 
neutrino, there is no interference from Eq.(A5) and thus no oscillation. The probability of 
remaining an electron neutrino remains one half without variation in the vacuum from the 
sun to the earth. In contrast, pure vacuum oscillations with maximal mixing (no MSW 
transition) leads to equal parts of each mass eigenstate; the neutrinos oscillate from pure 
electron neutrino to pure other flavor neutrino on the path from sun to earth. However 
the phase averaged probability of an electon neutrino reaching the earth is still one half. 



Guth, Randall, and Serna ||2^ have pointed out the relevance of this difference for matter 
oscillations in the earth: there can be a day-night effect, even for the case of maximal mixing 
if there has been an MSW transition in the sun. Appendix A comprises a short digression 
on this point. The treatment in Appendix A assumes phase averaging over the distances 
involved, due to the larger Am^ values that would come into play in a possible day night 
effect. Here we are interested in the phase of the vacuum oscillation, since that is our signal. 

The onset of MSW conversion in the sun with larger Am? can be investigated numerically 
by utilizing a piece of computer code adapted from a previous investigation [^]. The rate of 
^Be electron neutrinos emerging from the sun (again normalized to unity for no oscillations) 
takes the form 

RiArr?)=A + Scos- ^ , (7) 

^ ^ (0.862)0.00248' ^ ' 

with X the distance from the surface of the sun plus some constant. For maximal mix- 
ing: A = 0.5, and B is 0.5 for vacuum oscillations but B vanishes for complete adiabatic 
conversion. 

The top panel in Figure 2 shows how the magnitude of the oscillation for maximum 
mixing is reduced with increasing Am^ even though the phase averaged mixing remains 
a constant. The filled circles are the values of 2B calculated numerically for sin^ 2^^ = 1. 
The solid line through these circles approximates 2B by exp[— 1.583 (Am^)^]. Except for 
the small oscillations beyond Am^ = 10~*, this "MSW damping" is well represented by the 
gaussian factor. The solid line at .5 represents A for maximal mixing. The dotted lines are 
the corresponding quantities 2B and A for sin^ 29 = 0.9, just for comparison. 

Incorporating "MSW damping" Eq.(6) then becomes 



i?(0,Am2) = 1 - 0.395sin2 2^ (8) 
X 1 -exp(-1.583(Am2)2) |^c(- 



Am^ \ / 4396Am2 
cos 



+ .0167 COS 0/ VI + .0167 cos 

Am^ \ / 4396Am2 
+Is r sin 



,1 + .O167cos0y Vl + .0167cos( 
This is the expression that we use in the calculations to follow. 



Eq. (8) may also be written in the form 
i?(0,Am2) = 1 - 0.395sin2 2^ (9) 



X 



/ Am \ 

l-exp(-1.583(Am^)2) iJ - cos 

VI + .O167cos0/ 



1 + .0167 cos ( 



where 



S = arctan (10) 

-'c 

and 

lb = Ic sec 6. (11) 

For the purpose of illustration we ignore the cos dependence in the temperature damping 
factor and consider /;,(Am^). The bottom panel of Figure 2 shows /^(Am^) as the short- 
dashed line and repeats the gaussian MSW damping factor from the above panel as the 
long-dashed line. The solid line is the product of the two, the overall damping factor 
including temperature spreading and MSW damping. It is clear that there is a complete 
damping out of the oscillations at Am^ ~ 10^^ (sV)^, and that this broadening averages 
out the phase of the oscillations at higher values of Am^. 

Finally one should note that the source broadening is insignificant because of the follow- 
ing. It turns out that the SSM density [^] of ^Be neutrinos produced as a function of the 
solar radius is very close to a gaussian function of the sun's radius, 



Ws{r) = \l^e-^^\ (12) 



TT 



where f3 = 13.166 x 10^, and r is the distance from the center of the sun in units of the 



distance from the earth to the sun. However, as has been pointed out p7|JT8|] , the oscillations 
effectively start not at the source but at the level crossing point. For the present maximal 
mixing case the level crossing point is at surface of the sun. The neutrinos originating off 
the sun's axis in the direction of the earth will have a slightly larger distance to travel due 
to the curvature of the sun's surface. The gaussian density Eq.(12) leads, in a very good 
approximation, to a source spreading density of the form 



Ws{z) = fie 



(13) 



where /i = 1.23 x 10^ or twice (3 times the ratio of the sun's radius to the mean earth- 
sun distance, and z is the distance from the point on the sun's surface closest to the earth 
toward its center. If this small source broadening were the only cause of damping, then by 
an analytical treatment paralleling that leading to Eq.(9), one would find a source damping 
factor 



The long and short dashed line in the bottom panel of Fig. (2) represents /(Am^). Obviously 
/(Am^) only starts to deviate from unity at the rightmost part of the plot (at Am^ = 5). 
Thus, source broadening is insignificant for our region of interest. 

The region that we will investigate spans the range from Am^ ~ 10^^'', the "just so" 
region for ^B, up to Am^ ~ 10"^, where the broadening averages the phase. As noted above 
and in Appendix A, one might in principle begin to see a day-night effect with the onset 
of MSW damping. In fact there would be a sizable day-night effect for maximal mixing at 
Am2 ~ 10-^ [|9],|g (the so called "Low" MSW solution). 

Figure 3 shows R{4>, Am^) for maximal mixing sin^ 2^^ = 1 beginning at the low end 
with Am^ = 0.01 (again in units of 10^^(6^)^). Note from Eq.(8) that the overall phase of 
the cosine factor depends on Am^ and that this phase changes by 180° when Am^ changes 
by 7r/4396 = 0.000714. This phase sensitivity is illustrated Figure 3 where for each panel 
in addition to the curve for the labeled value of Am^ there are also curves for that value 
plus the appropriate increments to shift the overall phase by 90°, 180°, 270°, and 360°. 
Figure 4 shows the increasing frequency of the oscillations of R{(f), Am^) as a function of </> 
for the Am^ region of 10^^ to 10"*^. Note also the decreasing amplitude of the oscillations as 
they come close to being damped out by the temperature plus source broadening and MSW 
damping at 0.8 x 10^*^. 



/(Am^) 



1 



1 



(14) 
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III. OSCILLATIONS: ANALYSIS OF THE SIGNAL 



With such a rapid oscillation period throughout the year seen especially in Figure 4 (on 
the order of several months to several days), one anticipates that for such values of Am^ 
there would be insufficient statistics at an experiment like Borexino for a pattern to be 
obvious. However, in what follows we will investigate how a Fourier type analysis of data 
from such experiments could give evidence of a phased oscillation and thereby determine 
the value of Am^ if it lies in this range. Fourier analyses of ^Be solar neutrino data have 
been previously proposed [^,^], but what follows here is a somewhat different approach. 



Since Borexino is a detector rather than radiochemical experiment, it records the in- 
formation on when each count was recorded and thereby the distance of the detector to 
the sun L incorporated as (j) in Eq.(8). We suggest analyzing data from such experiments 
by effectively integrating data with a factor exp [i 4396Am^/(l + .0167 cos 0)] and varying 
Am^ over the range ~ 10^^'^ — 10^^ to look for a signal. 

To test whether the Am? can be determined by such a method, Monte Carlo data sets 
have been simulated in the following way. Random numbers are generated uniformly for 
from to 27r in order to cover the year and make use of Eq.(8). In order to weight events 
according to what would be expected from Eq.(8) with a specific Am^, a second random 
number between and 1 is then generated for each cj) and a count is generated if the random 
number is less than R{(f), Am'^). This is the data set: the collection of specific angles, {4>i}, 
at which single events are recorded durin a year. 

Figure 5 shows a sample analysis of a data set generated from 15000 Monte Carlo at- 
tempts for Am? = 0.3 in our units. The top panel shows the expected oscillation pattern, 
R{(f), 0.3). From Eq.(8) one would expect about 9075 data points to lie below the curve from 
15000 random attempts, and in fact a set of 8993 data points {(pi} were generated in this 
sample. The number of data points in this sample corresponds roughly to a year's running 
time at Borexino. For analysis one might first consider a Fourier type transformation on the 
set {(pi} 
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/ Am2, Am2 = - ^ cxp 

n^r^^ Vl + .O167cos0, 



(15) 



27ryo ^ ^ ^ Vl + .0167cos(/). 



The solid curve in the middle panel displays |/(0.3, Am^)| and it shows a discontinuity in 
pattern near Am^ = 0.3. If there were no phase oscillation then one would expect / to 
approach a function that we will call J(Am^) 

J(Am,2, = — / exp ^ (16) 

^ ^-Kh ^ Vl + .O167COS0; ^ ^ ' 

The dotted line in the middle panel displays |J(Am^)|. This suggests that we subtract 
off the Bessel function like behavior J(Am^) contained in /(Am^, Am^) (which tends to 
obscure the signal) , create a new function 

is:(Am^ Am^) = /(Am^ Am^) - J(Am^), (17) 

and use |i^(Am^, Am^)| to analyze our data set |i^(0.3, Am^)| is displayed in the 

bottom panel. The signal of Am^ = 0.3 is unambiguous. 

Figure 6 shows that an unabiguous signal would be obtained with about one year's 
Borexino statistics for Am^ in the range 0.1-0.5. For Am^ = 0.6 there is a signal at 
Am^ = 0.6 that might be a little ambiguous with only one year's statistics, but it retains 
its shape with increasing statistics; a secondary peak seen at about Am^ = 0.75 with one 
year's statistics goes away with the higher statistics. There is no signal apparent for the 
Am^ = 0.8 case, as one would expect from looking at the corresponding curve in Figure 4. 

Figure 7 shows the extraction of the signal for Am^ an order of magnitude lower. Curves 
correspond to values in Figure 3. Below Am^ = 0.04 this particular analysis starts to become 
ambiguous. However the slow rate of variation with ^ makes direct comparison with the 
patterns seen in Figure 3 practical. At Am^ = 0.01 the analysis is comphcated by the large 
change in magnitude of the rate throughout the year with a small increment in the value of 
Am^. 
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IV. DISCUSSION 



Based on the SNO result it now seem likely that the solar neutrino puzzle has been 
solved. It is not a deficiency in the standard solar model that is being observed but new 
physics. Electron neutrinos are oscillating into some combination of fi and r neutrinos. 
Exactly how this happens is perhaps not yet clear, whether by one of the MSW solutions 
or some vacuum mixing solution. In the previous sections of this paper it has been shown 
that if the solution to the solar neutrino puzzle happens to be maximal mixing in the mass 
range Am^ =~ 10^^*^ — 6 x 10~^(eV)^, then a proper analysis of a successful ^Be neutrino 
experiment should be able to unambiguously determine Am^. Not seeing a Am^ signal in 
this mass range would climate a region of Am^ for large mixing angle. 
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APPENDIX A: THE DAY-NIGHT EFFECT IN THE LIMIT OF MAXIMAL TWO 

NEUTRINO MIXING 

Guth, Randall, and Serna have pointed out that there can be a day-night effect, 
even for the case of maximal mixing. What follows is a compact explication of this point 
with emphasis on the limits of no MSW and maximal MSW effect in the sun. 

The general form for two mass eigenstates in two neutrino mixing is 

|z/i>= cos9\h'e> + sm9\h'x> (Al) 

and 

|z/2>= —sm9\h'e> + cos9\h'x> , (A2) 
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where \vx> is presumed to be some linear combination of \i'pi> and |i^r>. Conversely 

\Vf. >— cos^li^i > — sm.9\v2 > (A3) 

and 

\vx>— sni9\vi> + cos9\v2> . (A4) 

In free space mass eigenstates propagate independently. A mixed mass state \i'{t) > then 
has the form 

\p{t)>= e-''"?*/'^Ai|z/i> + e-^"*'*/'^A2|t/2> . (A5) 

The probability Pg that a neutrino born in the sun is an electron neutrino when it reaches 
the earth is then 

Ps = cos^^Pi + sin^^Pa (A6) 

with Pi 2 = 1^1,2 P the average probability of a mass one or mass two eigenstate arriving at 
the earth where the phase has been averaged by the distances involved. Since Pi = 1 — P2 
this may also be written equivalently 

^ (COS^^ - Ps) 



(AT) 



The probability P that an electron neutrino born in the sun will be an electron neutino 
after passing through the sun, traveling to the earth, and then passing through the earth is 
simply 

P = P2P2e + (1 - P2)(l - P-Ze) (A8) 

where P2e is the probability that a mass 2 neutrino entering the earth emerges at the detector 
as an electron neutrino. Making use of Eq. (A7) this becomes 
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This expression is the Mikeyev-Smirnov expression for the day night effect, trivially 



transformed [^] to be most transparent in various limits. 
The maximum value for Pg occurs for vacuum oscillations 

1 



Ps = l- ^sm^2e (AlO) 



and 

fl-2R 



cos2^ 



-cos20. (All) 



The minimum value for Ps occurs for complete adiabatic MSW conversion to a pure mass 
eigenstate |z/2>. In this case from Eq.(A6) 

Ps = sin^e (A12) 

and 

^ TT^ = 1- (A13) 

cos2^ ^ ^ 

Thus as sin 26 goes to 1 (maximal mixing) there is no day-night effect for vacuum oscil- 
lations and a maximum effect possible in the case of complete adiabatic conversion. 
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FIGURES 





0.01 0.1 1 

2 -8 2 

Am in units of 10 (eV) 

FIG. 2. Top panel: Constant A and oscillating B part of electron neutrino rate emerging from 
the sun for sin^ 2^ = 1 (solid line) and for sin^ 2^ = 0.9, (dotted line); Bottom panel: Total 
oscillation damping factor (solid line) and partial damping factors (see text). 
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90 180 270 360 90 180 270 360 



FIG. 3. Rate of Be neutrinos detected by electron scattering for maximal mixing. The number 
on each panel gives Am^ for the solid line. The long and short dashed line is Am^ plus approx- 
imately 0.00035, the short dashed line Am^ plus 0.0007, the long dashed line Am^ plus 0.00105 
and the dotted line Am^ plus 0.0014. 
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90 180 270 360 90 180 270 360 

FIG. 4. Rate of ^Be neutrinos detected by electron scattering for maximal mixing. The number 
on each panel gives Am^ for the solid line. 
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90 180 270 360 




0.5 1 

2 -8 2 

Am^ in units of 10 (eV) 

FIG. 5. The top panel gives the rate of ^Be neutrinos coming from the sun for maximal mixing 
and Am^ = 0.3. The solid line in the middle panel gives the absolute value of the Fourier analysis 
of the distribution |7(0.5, Am^)| and the dashed line the absolute value of the analysis of the 
constant average of the distribution | J(Am^)|. The bottom panel shows the absolute value of the 
difference of the Fourier analyses of the distribution and the constant average of the distribution 
|K(0.5, Am^)|. This is the quantity that extracts the Am^ signal. 
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0.5 1 0.5 1 



Am^^ in units of 10 \eV)^ Am^^ in units of 10 \eV)^ 

FIG. 6. Extraction of the the signal |i^(Am^, Am^)|. The number on each panel gives the 
value of Am^. The dotted line lines represents 15000 tries or about 9000 events. The solid lines 
represent four times the statistics: 60000 tries or about 36000 events. 
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Am^ in units of 10 (eV) Am^ in units of 10 (eV) 

FIG. 7. \K {Am^ , Aml)\ as in Figure 6, but for lower values of Am^. The location of the peak 
becomes less well defined with decreasing Am^. As in Figure 3 the long and short dashed line 
is AmP' plus approximately 0.00035, the short dashed line Am^ plus 0.0007, the long dashed line 
Am? plus 0.00105 and the dotted line Am^ plus 0.0014. 
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